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Abstract. To identify the nature of the newly observed charged resonance Z c ( 4200) + , we study its 
hadronic decays Z c ( 4200) + — > J/ipn + , Z c (4200) + — >• r/ c p + and Z c (4200) + —> D + D*° as a charmonium-like 
tetraquark state. In the framework of the QCD sum rules, we calculate the three-point functions and ex¬ 
tract the coupling constants and decay widths for these interaction vertices. Including all these channels, 
the full decay width of the Z c (4200) + state is consistent with the experimental value reported by the Belle 
Collaboration, supporting the tetraquark interpretation of this state. 

PACS. 12.38.Lg Other nonperturbative calculations - 11.40.-q Currents and their properties - 12.39.Mk 
Glueball and nonstandard multi-quark/gluon states 


1 Introduction like states Zb(10610) and Zb( 10650) in the 7r ± T’(nS') and 

hb^ mass spectra in the T(5S) decay [TT] , 

Recently, a new charged charmoniumlike resonance Z c (4200) + These newly observed charged states_have the exotic 
was observed by the Belle Collaboration [T] . It was ob- flavor contents ccud for Z c states and bbud for Zb states. It 


served in the Z c ( 4200) + —> J/ilm + process with the mass 
and decay width M = 4196^29^13 MeV and T = 370+7oii3 2 
MeV, with a significance of 6.2cr. Its preferred assignment 
of the quantum numbers is J p = 1 + . The G-parity of 
Z c (4200) + is positive. Thus, the quantum numbers of its 
neutral partner is I G J PC = l+l 4 . 

The family of the charged charmoniumlike states have 
become more abundant after the discovery of Z c ( 4200) + [X] 
and Z c (4050) [2]. Before this, the first member Z(4430) + 
was observed in the ifj{2S)'K + invariant mass spectrum in 
the process B° —> %Ij(2S)tt + K~ by the Belle Collabora¬ 
tion [3] and confirmed recently by the LHCb Collabora¬ 
tion g]. Later, Belle also reported a broad doubly peaked 
structure in the 7r + x c i invariant mass distribution, which 
are called Z(4050)+ and Z{ 4250)+ [5]. Several other simi¬ 
lar charged states were observed in last two years. In 2013, 
the BESIII Collaboration reported Z c (3900) + in J/i/jtt + fi- 


is natural to understand them as different manifestations 
of four-quark states: hadron molecules mmmmmrn 
CEB], tetraquark states mmm, or many other configu¬ 
rations [2211231123] . For example, Z{ 4430) + was described 
as a D*D\ molecular state in Refs. [25ll26ll27ll28] and a 
tetraquark state in Refs. [231130113lj : the Z c (3900)+ was 
speculated to be a molecular state in Refs. [3211331133] : the 
Z c ( 4025) + was interpreted as a D*D* molecular state in 
Refs. [3S]; the Z&(10610) and Zb( 10650) were studied as 
BB* and B*B* molecular states in Refs. [33]. One can 
consult Refs. [37113311331133113X1132] and references therein 
for recent reviews of these charged resonances. 

Being composed of a diquark and antidiquark pair, a 
hidden-charm tetraquark state can decay very easily into 
a pair of open-charm D mesons or one charmonium state 
plus a light meson through quark rearrangement, imply¬ 
ing that tetraquark states should be very broad resonances 


nal states in the process Y (4260) 


J/lp7T + TT 


was also observed by Belle [7] and confirmed in CLEO 
data [5j. The BESIII Collaboration also observed Z c (4025) ± 
in the recoil mass spectrum in the e + e _ —> (D*D*)^ 1 
process [3] and Z c (4020) ± in the h c n ± mass spectrum in 
the process e + e _ —> h c n + n~ - m ■ Moreover, the Belle 
Collaboration also observed two charged bottomonium- 
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[6]. Z c (3900) + while the experimental XYZ states are usually quite nar¬ 
row, such as Z c (3900) + [lOEM] an d ■Z’ c (4025) + [311X3]. How¬ 
ever, the experimental width value of the Z c (4200) + jXj 
is broad enough to be a good tetraquark candidate. In 
Ref. [33], ir c (4200) + was studied as a tetraquark state by 
considering the color-magnetic interaction. In Ref. [44] . 
the authors tried to search for Z+ exotic states in lattice 
QCD. However, they found no convincing signal for Z+ 
state below 4.2 GeV. 
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The hidden-charm tetraquark states with J PC = l -1 
has been studied using the method of QCD sum rule 
in Refs. [MirniTOirnpTRipnTifTiTii . We have also done sim¬ 
ilar QCD sum rule studies in Refs. EEE2, in which the 
extracted mass was found to be consistent with the ex¬ 
perimental value of the Z c (4200) + mass. In this work, 
we will study the hadronic decays of the if c (4200) + as a 
tetraquark state in QCD sum rules. The three-point func¬ 
tions for the ZcJ/ipTr, Z c r] c p and Z C DD* vertices will be 
studied to calculate the corresponding coupling constants 
needed to extract the decay widths. 

This paper is organized as follows. In Sec. El we study 
the three-point functions for the Z C J/if) 7r, Z c r] c p and Z C DD* 
vertices. We will calculate the OPE series up to dimension 
five condensates. Then we compute the coupling constants 
and the decay widths for these channels. Finally, we give 
a short summary and discuss the possibility of searching 
for such exotic resonances decaying into q c charmonium. 

2 QCD sum rules and three-point correlation 
function 

In the past several decades, QCD sum rule has proven to 
be a very powerful non-perturbative approach to study 
hadron properties such as masses, magnetic moments and 
coupling constants, associated with the low-lying baryons 
and mesons [53l(54l(55ll561l57l . Recently, this method was 
used to yield predictions on the spectroscopy of the new 
hadron XYZ states 

To calculate the decay width of the Z c ( 4200) meson 
into two hadrons, one needs first to study the three-body 
coupling vertices Z C (4200)AR, where A, B denote the de¬ 
cay products. In QCD sum rules, we consider the three- 
point correlation function 

n^(P:P',q) (1) 

= J d A xd A y e ip '' x e iq ' v (0\T[J A (x)J B (y)J^(0)}\0), 

where J l f c is the interpolating current for the Z c ( 4200) + 
meson while J A and J B are the currents for the final 
states A and B, respectively. In this paper, we consider the 
Z’ c (4200) + meson as a charmonium-like tetraquark state. 
The corresponding tetraquark current is given by 

J u c = u T a C l5 c b (d alv Ccl + d b ^Ccl) (2) 

-u^C^ v c b {d a ^Ccl +d b j 5 Ccl), 

in which the subscripts a, b are color indices, and u , d and 
c represent up, down and charm quarks, respectively. C 
is the charge-conjugation matrix. We have studied this 
charmoniun-like tetraquark scenario and the extracted mass 
is around 4.16 GeV m consistent with the observed mass 
of the Z c (4200) + meson [I]. This current can couple to the 
Z c (4200) + meson via 

(0\J?\Z c (p)) = f z e„(p), (3) 

in which e w {p) is a polarization vector and fz is the cou¬ 
pling constant of the current to the physical state. 


The Z c (4200) meson can decay into several different 
channels such as hidden-charm decay mo_des J/ipn + ,r] c p + 
and open-charm decay modes D + D*°, D°D* + . Such de¬ 
cay properties are similar to those for the charmonium- 
like state Z c (3900). Assuming Z c (3900) to be a tetraquark 
state with the same quantum numbers as the Z c (4200), 
the hadronic decay modes of Z c (3900) to J/'i/jn + , q c p + , D + D*° 
and D°D* + were studied in Ref. m ■ Building upon these 
methods, we will study the same decay channels for the 
Z c (4200) tetraquark state to estimate its decay width. 

2.1 Decay mode Z+(4200) —> J/i/jit + 

In this subsection, we study the hidden-charm decay Z p (4200) 
J/ijm + , in which mode the Z p (4200) meson was observed. 

To calculate the three-point function in Eq. © for the ver¬ 
tex Z+(4200) J/^7r + , we need the interpolating currents 
for J/ip and 7r mesons 

Jt = Cal^C a , (4) 

J W = da”f5U a , ( 5 ) 

which can couple to the J/ip and n mesons respectively 
via the following relations 

(o\j*m(j/)) = ( 6 ) 

<01 J>(g)) = /; = ^M (7) 

J 7T 

where and /(. are coupling constants. Using these re¬ 
lations and Eq. ©, we can write down the three-point 
function in the phenomenological side 

n^(p,p',q) (8) 

= J d 4 ‘xd i y e ip ' x e iq v <0|T[ (x)J v {y)J^ (0)]|0) 

( P'A \ ( PvPv'\ .. 

- gz^'u' to 2 J [9^' m | ) x 

__ + 

( p 2 — m| + ie)(p' 2 — m ^ + ie){q 2 — m 2 + ie) 

= _ gz^g^hi-if'^fz _ x 

(p 2 — m 2 z + ie)(p' 2 — m ^ + ie)(q 2 — m 2 + ie) 

q^p’u + g^9v _ PyPv_ p' ■ gjp'yP'v + \ 

m| rri^ rn z rn '^ J 

+..., 

in which represents the contributions of all higher 
excited states. We have used the relation p = p' + q in the 
last step. The coupling constant (form factor) gziim(q 2 ) is 
defined via 

£ = gz^g^Z+TT-^p,, + h.c.. (9) 

In Eq. ©, the three-point function n^(p,p\q) is diver¬ 
gent at q 2 = 0 when we take the limit m, r = 0. To simplify 
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the calculation in OPE side, we establish a sum rule at the 
massless pion-pole, which was first suggested in Ref. 54] 
for the pion nucleon coupling constant. 

At the quark-gluon level, the three-point function in 
Eq. (JTJ) can be evaluated via the operator product expan¬ 
sion (OPE) method. We insert the three currents, Eqs. ©, 
Q and ([5]), into the three-point function, Eq. ©, and do 
the Wick contraction: 


(0\T[J*(x)J”(y)JZH0)]\0) (10) 

= -T^iSab’ixhsS'ziyhsS'Jbi-yhuS^^-x)^,,} 
-^[Sab'i^Sba'iyhsS^bi-y^S^^-x)^] 

-^[Sab'ixhvSC’iyhbSvbi-yhsSa'ai-xhv.}, 


where the subscripts a, 6, a' and b' are color indices, and 
the superscripts u, d and c denote the quark propagators 
for up, down and charm quark, respectively. Throughout 
our evaluation, we use the coordinate-space expression for 
the light quark propagator and momentum-space expres¬ 
sion for the heavy quark propagator [541(59] : 


iS ab( x ) = ^~rX 


7 \ n 1 

- S -§m + ^(MoGq), 


xcr 


fiv\ 


iS c ab {p) = 


iSab(j/> + rn c ) 

p 2 — m 2 


i \ n 1 

1 A ab rrn 1 

2 ^(p 2 -m 2 ) 2 


(j> ± m c ) + (j> + m c ) 


( 11 ) 


( 12 ) 


where m c is the mass of the charm quark. We neglect the 
chirally-suppressed contributions from the current quark 
masses ( m q = 0 in the chiral limit) because they are nu¬ 
merically insignificant. In Eq. m, the light quark prop¬ 
agator is defined as iS^ b (x) = C(iS^ b ) T C in which T rep¬ 
resents only the transpose operation to the Dirac indices. 
As indicated above, we will pick out the 1 /q 2 terms in the 
OPE series and work at the limit q 2 —> 0. We note that 
this is the assumption used in Ref. |54j . and then we can 
establish a sum rule by comparing with the three-point 
function expression Eq. © at the hadron level. 

On the phenomenological side in Eq. ©, there are five 
different tensor structures g^, q^p' u , q^p'^ and p'^p'v 
On the QCD side, we evaluate the three-point function 
and spectral density up to the dimension five terms. In 
addition to the perturbative term, we calculate the quark 
condensate, the gluon condensate and the quark-gluon 
mixed condensate for the power corrections. The Feyn¬ 
man diagrams for these terms are shown in Fig. [H In our 
result for the OPE series, the g structure contributes to 
all expansion terms including the perturbative part, quark 
condensate, gluon condensate and quark-gluon mixed con¬ 
densate. Other tensor structures contribute just some of 
these terms in the OPE at leading order. For example, the 
q^qv and q^p'^ structures appear only in the gluon conden¬ 
sate while p'nPv appears in the perturbative term and the 
gluon condensate. The structure q^p'^ gives no contribu¬ 
tions to the perturbative term. 


To obtain the greatest number of terms in the OPE se¬ 
ries, we therefore study the gy v structure in the following 
analysis. Finally, we obtain the spectral density propor¬ 
tional to 1 /q 2 in the g^ structure 


(g 2 GG)(s + 2m 2 c ) / 4mf 
192 tt 4 V s 


(13) 


where we find that only the gluon condensate gives con¬ 
tributions to the spectral density at order 1/q 2 . There are 
three kinds Feynman diagrams for the gluon condensate 
in Fig. [T] in which the third one is color-connected and the 
former two are color-disconnected m- In our caculation, 
the spectral density in Eq. m comes from only the first 
color-disconnected diagram in Fig. [T| The second color- 
disconnected diagram and the color-connected diagrams 
give no contribution to p(s). 

To establish a sum rule for the coupling constant gzi/jv, 
we assume p 2 = p' 2 = P 2 in Eq. © and then perform the 
Borel transform (P 2 —> M |) to suppress the higher state 
contributions. For the g^ structure, we arrive at the sum 
rule 


gZiPTr{SQ, M|)|q2^q 
1 


mP — mf, 


U(-if*)fz e ~ m * /M B - e- m z/ M B 


(14) 

r s ° , 9 

/ p(s)e- s/M zds, 
J 4 m? 


in which Q 2 = — q 2 and so is the continuum threshold 
parameter for the Z c (4200) meson. 

To perform the QCD sum rule numerical analysis, we 
use the following values of quark masses and various con¬ 
densates [MMimi]: 


m c (p = m c ) = m c = (1.275 ± 0.025) GeV, 

(qq) = -(0.23 ± 0.03) 3 GeV 3 , 

{qg s a ■ Gq) = -M${qq) , (15) 

Ml = (0.8 ± 0.2) GeV 2 , 

{g 2 s GG) = (0.48 ± 0.14)GeV 4 , 


in which the charm quark mass is the running mass in 
the MS scheme. Note that there is a minus sign implic¬ 
itly included in the definition of the coupling constant g s 
in this work. We use the following values of the hadron 

rriqp = (3.097 ±0.011) GeV, /p = (1.288 ± 0.037) GeV 2 , 
m n = 139.6 MeV, U = 133 MeV, (16) 

m Zc = (4.196ig;gg) GeV, f Ze = (6.9 ±0.4) x 10“ 3 GeV 5 , 

in which the coupling parameter fz c is determined by the 

mass sum rules in Ref. 1511- 

In Eq. m, the coupling constant gz^{q 2 ) depends on 
the continuum threshold value so and the Borel mass Mb 
in the limit q 2 —> 0. We use the continuum threshold value 
so = 21 GeV 2 , which is the same value as used in the mass 
sum rule in Ref. l5ll . Using this value of so and the param¬ 
eters given in Eqs. m and m, we show the variation of 
the coupling constant gz^-K with the Borel parameter Mg 
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Fig. 1. Feynman diagrams for the leading order contributions of the three-point function. Wave lines depict vector and axial- 
vector mesons, dashed lines pion, solid lines quarks and curly lines gluons. Graphs related by symmetry are not shown. 



Ml [GeV 2 ] 

Fig. 2. Coupling constant gz^i r dependence on the Borel pa¬ 
rameter Mg. 


where gABC is the coupling of the three-point vertex ABC 
and p*(mA,mB,mc ) is defined as 

^~2aAlA~~~2c?a? 

P {a,b,c) = ---(18 

2 a 

For the decay mode i?+(4200) -A J/ijjn + , we then can 
calculate its decay width using Eq. m 

r(Z+( 4200) -a J/^7T+) = (87.3 ± 47.1) MeV, (19) 

in which the dominate error source is the uncertainty of 
the gluon condensate. 


2.2 Decay mode Z+(4200) — y g c p + 


in Fig. [2j We find that the sum rule gives a minimum (sta¬ 
ble) value of the coupling constant gz^v = (6.27 ± 1.93) 
GeV with M| ~ 1.9 GeV 2 . The errors come from the un¬ 
certainties of the charm quark mass, the gluon condensate, 
hadron masses and hadron couplings. To calculate the de¬ 
cay width of a process A -A BC , we use the following 
relation given by Refs. mm-- 


We study the decay mode Z+ (4200) -A rj c p + in this sub¬ 
section. To calculate the corresponding three-point func¬ 
tion, we consider the following interpolating currents for 
r) c and p mesons 

J Vc = c a 7 5 c Q , (20) 

= dal^Ua, ( 21 ) 

with the current-meson coupling relation 


r(A -A BC) 

_ p*(m A , m Bl m c ) 

87 Tm 2 A 


9abc 


(17) 

/ 3 p*{m Al 'mB,rnc) 2 \ 

V rn? B ) ’ 


f TO 2 

(0\J^\Vc(p')) = fi (22) 

(0\Wq)) = M<l), (23) 


3 
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where f Vc and f p are coupling constants for r\ c and p, 
respectively. Then we can obtain the three-point function 
at the hadron level 

n™(p,p',q) ( 24 ) 

= J d i xd 4 ‘y e ip '' x e ig ' y (0\T[J^(x) J p (y) J*+(0)]|0> 

= _ 9Zr, cP (q 2 )fp(-if^ c )fz _ x 

(p 2 — m 2 z + ie)(p' 2 — m 2 c + ie)(q 2 — m 2 + ie) 

qyp'p + p'pP'v _ qp,gv p' ■ q(qpqv + qp.p' v ) \ 
m 2 z m 2 m 2 z m 2 J 

+•••> 



where the coupling constant gzri c p(q 2 ) is defined via 

£ = gzr, cP g^ Z+r] c p~ + h.c .. (25) 

At the quark-gluon level, we calculate the three-point 
correlation function 7J^ p (p,p', q) by considering diagrams 
similar to Fig.[Q As outlined in Ref. m , the coupling con¬ 
stant varies slowly with the Euclidean momentum Q 2 = 
—g 2 , and hence for sufficiently large Q 2 it is only necessary 
to extract the 1 /Q 2 term in Eq. (1241) for the three-point 
correlation function njfc p (p,p',q). Thus, for an appropri¬ 
ate range of Q 2 we keep the invariant function propor¬ 
tional to l/Q 2 in the g structure. In this assumption 
m, we obtain the spectral density on the OPE side 

m c (qg s a ■ Gq) {g 2 GG)s\ f 4m 2 "^ 
-(hr 2 - + ^92^J 



in which only the gluon condensate and quark-gluon mixed 
condensate give contributions to the spectral density in 
this order. Similarly to the spectral density in Eq. C3D, 
the color-connected diagrams in Fig. |T] give no contribu¬ 
tions to this spectral density in Eq. (1261) . Then the sum 
rule for the coupling constant gz Pc p can be established by 
assuming p 2 = p' 2 = P 2 and performing the Borel trans¬ 
form (P 2 -)> M%) 



in which Q 2 = — q 2 . 

To perform the numerical analysis, we use the following 
parameters for the p c and p mesons mm 


m Vc = (2.980 ± 0.001) GeV, f Vc = 0.35 GeV, 
to p = (775.26 ± 0.25) MeV, f p = 157 MeV. (28) 

In Fig. [3j we show the variation of gzr/ cP (Q 2 ) with 
the Borel parameter for Q 2 = 8 GeV 2 . We choose m 2 <C 
Q 2 ~ m 2 c so that the m p can be ignored and the OPE 



Fig. 3. Variation of the coupling constant gzri cP (Q 2 ) with the 
Borel parameter Ml for Q 2 = 8 GeV 2 . 


is valid. It shows that the minimum value of gz-q a p ap¬ 
pears around Ml = 1.6 GeV 2 . However, the coupling con¬ 
stant is required at the pole Q 2 = — m 2 , where the QCD 
sum rule is not valid. To obtain the coupling constant 
gz VcP (Q 2 = —m 2 ), we need to extrapolate the coupling 
constant gzr] c p{Q 2 ) from the QCD sum-rule region to the 
physical pole Q 2 = — to 2 . Following Ref. m > we use the 
exponential model to achieve this extrapolation 

gzp c p(Q 2 ) = gie~ 92Q2 , (29) 

where gi and 52 can be determined by fitting the QCD 
sum-rule result of gzri cP {Q 2 ), using Eq. (l29l) . In Fig. [4] we 
show the QCD sum-rule result of the coupling constant 
gz VcP (Q 2 ) with s 0 = 21 GeV 2 and Ml = 1.6 GeV 2 (dot¬ 
ted line). We fit this result by using the model Eq. (EH1) 
with the parameters g\ = 11.65 GeV and 52 = 9.93 x 
10 -3 GeV 2 . Then we extrapolate the coupling constant 
gz VcP (Q 2 ) to the physical pole Q 2 = —to 2 to obtain the 
coupling constant 

gzp cP (Q 2 = -m 2 ) = (11.72 ± 2.10) GeV, (30) 

in which the uncertainties of the charm quark mass, the 
QCD condensates, the hadron masses and the hadron cou¬ 
plings are considered to give the error of this prediction. 
Inserting this value into Eq. m, we can calculate the 
decay width of the process Z+ (4200) —> r] c p + 

r{Z+{ 4200) rj c p + ) = (334.4 ± 119.8) MeV. (31) 


2.3 Decay mode Z+(4200) -> D+D*° 

In this subsection,_we study the open-charm decay mode 
21+(4200) —> D + D*° with the interpolating currents for 


the D + and D*° mesons 

J — da^i^Ca’i 

(32) 

t D * 0 _ - 
J H — c a r ln u ai 

(33) 

with the current-meson coupling relations 

(0\J d+ \D+(p')) =f Dl 

(34) 

<0| jfV°(g)}=/W<z)> 

(35) 
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12 . 0 - 



100 m | 0.5 1.5 2.5 3.5 4.5 5.5 6.5 7.5 8.5 


Q 2 [GeV 2 ] 

Fig. 4. Variation of the coupling constant gz PcP {Q 2 ) with Q 2 . 
Data points shows the QCD sum-rule result with so = 21 GeV 2 
and Mg = 1.6 GeV 2 . Solid line gives the fit of the QCD sum- 
rule result through Eq. (|29]i and the extrapolation of the cou¬ 
pling constant gzrj cP {Q 2 ) to the physical pole Q 2 = —m 2 . 


where fn and fn* are the current-meson coupling con¬ 
stants for D and D* , respectively. The three-point func¬ 
tion can be written at the hadron level 

nj?„ D * (p,p',q) ( 36 ) 

= J d A xd A ye'v'-*e^y (0\T[J D (x)J°'(y)J?\0)}\0) 

_ _ gzpp- ( q 2 )fp * {—ifp)fz _ x 

(p 2 — m 2 z ± ie)(p' 2 — m 2 D ± ie)(q 2 — m 2 D , ± ie) 

QvP’n + P'^Pl q P q P p' ■ q{q,j.qv + q P p' p ) \ 

o 2 ' 2 2 I 

m z m D * m z m D * J 

+•••> 



where the coupling constant gzDD*{q 2 ) is defined via 

C = gzDD'g^Z+D-D; 0 ± h.c .. (37) 

However, the evaluation of the three-point correlation 
function 77 (p, p', q) is a bit different from the other two 
decay modes discussed previously. Considering diagrams 
similar to Fig. |T| we calculate the OPE series in the mo¬ 
mentum spaces. As discussed above, we keep the invariant 
function proportional to 1/ Q 2 in the structure at the 
OPE side 


,<gGg>( s ) = rn c {qg s a ■ Gq) 


6 - 


5m 2 


19 m 2 
s — m 


192tt 2 
(g 2 s GG) 

61447r 4 s 2 . 

— (15m 2 s 2 — 34m 4 s — to®) log ( 1 + 


(4s d + 19m 2 s 2 — 34m 4 s — 5 m®) 

2 ' 


m c 

Q 2 


(38) 


in which Q 2 = —q 2 . Here the quark-gluon mixed con- 
densated (qg s c • Gq) gives dominant contribution to the 
spectral density shown above. We find that for ( qg s a-Gq ), 
only the color-connected diagrams contribute to the spec¬ 
tral density shown above, which is different from the situa¬ 
tions in J/tjjTT and g c p channels. For the gluon condensate 



Fig. 5. Variation of the coupling constant gzDD*(Q 2 ) with 
the Borel parameter Mg for Q 2 = 20 GeV 2 . 



Fig. 6. Variation of the coupling constant gzDD*(Q 2 ) with 
Q 2 . Data points shows the QCD sum-rule result with so = 
21 GeV 2 and Mg = 3.4 GeV 2 . Solid line gives the fit of the 
QCD sum-rule result through Eq. (1291) and the extrapolation 
of the coupling constant gzDD * ( Q 2 ) to the physical pole Q 2 = 
-m 2 D ,. 


(g 2 GG), both the color-connected and color-disconnected 
diagrams give contributions. Assuming p 2 = p' 2 = P 2 
and performing the Borel transform (P 2 —> Mg) of the 
three-point function, the sum rule for the coupling con¬ 
stant gzDD * can be obtained as 


gzDD*(s 0 , Mg, Q 2 ) 


1 m 2 z — m 2 D 

fD* ( ifD)fZ e _m n/ M B — e - m z/ M B 



(39) 


We adopt the hadron parameters of D + and D*° from 
Refs. |37U47] 


m D + = (1869.61 ± 0.10) MeV, f D = (0.18 ± 0.02) GeV, 
m D ,o = (2006.96 ±0.10) MeV, / D * = (0.24 ±0.02) GeV. 

(40) 

In Fig. 0 we show the variation of gzDD*(Q 2 ) with 
the Borel parameter while Q 2 = 20 GeV 2 . In this situa¬ 
tion, the coupling constant increases monotonically with 
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the momentum Q 2 , and hence the result obtained below 
should be considered as an upper bound on the compar¬ 
atively small decay width in this channel. To perform 
the QCD sum-rule analysis, we adopt the Borel window 
3.0 < Mj, < 3.4 GeV 2 used in Ref. [511 f° r the mass sum 
rules of the same current. We do the fitting at M| = 3.4 
GeV 2 and So = 21 GeV 2 in Fig. [6] using Eq. (l29l) with 
gi = -1.73 GeV and g 2 = 1.77 x 1CT 2 GeV 2 . With 
these parameters, the model Eq. C9l) can fit the QCD 
sum-rule result very well. To obtain the coupling con¬ 
stant, we extrapolate the coupling constant gzDD*(Q 2 ) 
to the physical pole Q 2 = — m 2 D „. Considering the un¬ 
certainties of the input parameters, the coupling constant 
gzDD*(Q 2 = - m 2 D ,) is 

gzDD *(Q 2 = —m 2 D .) = -(1.86± 1.09) GeV. (41) 

Finally, the decay width of this decay mode can be evalu¬ 
ated via Eq. ED 

r(Z+( 4200) D + D*°) = (6.6 ± 6.4) MeV. (42) 

In addition, the open-charm decay mode Z+ (4200) —> 
D°D* + should also be studied. Considering the SU(2) 
symmetry, the three-point function for the vertex 
Z+ (4200 )D°D* + is exactly the same as that for the ver¬ 
tex Z+ (4200) D + D*°. Using the hadron parameters of D° 

and D*+ [37]: 

m D o = (1864.84 ±0.07) MeV, 
m D *+ = (2010.26±0.07) MeV, 

we obtain the same values of the coupling constant and 
the decay width for Z+(4200) -± D°D* + as those for 
Z+ (4200) -> D+D*°. 

3 SUMMARY 

In summary, we have studied the three-point functions of 
the processes Z c (4200)+ —> J/^7r + , Z c (4200) + —> rj c p + 
and Z c ( 4200) + —> D + D*°, considering Z c (4200) + as a 
hidden-charm tetraquark state. We calculate the three- 
point functions by including the perturbative term, quark 
condensate, gluon condensate and quark-gluon mixed con¬ 
densate. 

To perform the QCD sum rule analysis, we expand the 
three-point functions in QCD side with respect to Q 2 and 
isolate the terms proportional to 1/Q 2 from the g^ ten¬ 
sor structure. Only the gluon condensate and the mixed 
condensate give contributions to the three-point functions 
after this procedure. This approach is a modification of 
the QCD sum-rule analysis of the Z c (3900) + decay width 
in Ref. [47] where the three-point functions were evalu¬ 
ated by considering the color connected (CC) diagrams 
associated with a different tensor structure resulting in 
only mixed condensate contributions. Our predictions of 
the decay widths are 

r(Z c ( 4200)+ -± J/ipn + ) = (87.3±47.1) MeV, 
r(Z c ( 4200)+ r/ c p+) = (334.4 ± 119.8) MeV, 

F(Z C (4200)+ -± D + D*°) = (6.6 ± 6.4) MeV, 
r(Z c { 4200)+ D°D* + ) = (6.6 ± 6.4) MeV. 


Thus, the full decay width for Z c { 4200) + is predicted as 

r Zc{ 4200 )+ = (435 ± 180) MeV, (43) 

which is in agreement with the experimental value of Z' c (4200) + 
width from the Belle Collaboration pQ. It is found that 
the branching fraction into J/tjjiT channel is about 24.7%, 
which is slightly suppressed compared to 71.6% for the 
r] c p channel. The study of Z c —> rj c p decay may pro¬ 
vide useful insights on the nature of the newly observed 
charged Z c states, helping to discriminate the molecule 
and tetraquark interpretations of the charged state fam¬ 
ily [53]. 

The study of the three-point function sum rules gives 
support to the tetraquark interpretation of the newly ob¬ 
served Z c ( 4200) + state. This conclusion is consistent with 
the result obtained from the mass sum rules in Ref. [51] . 

The branching ratio predictions of J/i/jtt, r] c p, D + D*° and 
D°D* + channels will be helpful for future experimental 
studies. 
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